Infrared behavior and gauge artifacts in de Sitter spacetime: The photon field 
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We study the infrared (long-distance) behavior of the free photon field in de Sitter spacetime. 
Using a two-parameter family of gauge-fixing terms, we show that the infrared (IR) behavior of the 
two-point function is highly gauge-dependent and ranges from vanishing to growing. This situation 
is in disagreement with its counterpart in flat spacetime, where the two-point function vanishes in 
the IR region for any choice of the gauge-fixing parameters. A criterion to isolate the "physical" 
part of the two-point function is given and is shown to lead to a well-behaved two-point function in 
the IR region. 
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Introduction. The study of quantum field theory in 
the de Sitter background is of paramount importance to 
the understanding of the early Universe as well as its 
present accelerated expansion. One of the most strik- 
ing, yet poorly understood, aspects of de Sitter quantum 
field theory (QFT) is the so-called infrared problem. The 
simplest and most serious manifestation of this problem 
being the nonvanishing of the correlation functions for 
largely - spacelike and timelike -separated points. This 
is for instance a well-known fact for the massless mini- 
mally coupled (mmc) scalar and the graviton fields. In 
this Letter we study the IR behavior of the photon field 
in de Sitter space and show that it exhibits similar IR 
pathologies. More importantly we are able to show that 
these are purely gauge artifacts. 

The organization of the Letter is as follows. First we 
compute the correlation function using the R\ gauges 
A(V a A a ) 2 /2. Then we generalize the computation to 
the larger family of higher derivative gauge-fixing terms 
\{V a A a ) 2 /2 + £,R 2 {V a A a )U{V a A a )/2. Finally we pro- 
pose a decomposition of the two-point function into a 
physical and a gauge part and we show that the physical 
part is vanishing in the IR region. 

The d-dimensional de Sitter spacetime Xd of "radius" 
R can be identified with the real one-sheeted hyperboloid 
in the d+1 Minkowski spacetime Md+i 



X,, 



{x G 



V 3 



R 2 } 



where r/^ is the "mostly plus" flat metric 
diag(— 1, +1, • ■ ■ ,+1). Let fjt(x,x') denote the geodesic 
distance between the points x and x' and g ab be the 
d-dimensional de Sitter (or d-dimensional Minkowski 
if R —> oo) metric. We will compute the Wightman 
two-point function W aa >(x,x') = (rt\A(x)A(x')\n) in 
different gauges. In the previous formula \Q) designates 
the Euclidean (Bunch-Davies) vacuum state. This 
vacuum being de Sitter invariant, W aa >(x,x') is a 
maximally symmetric bi-tensor. As shown in detail in 
[H, all such bi-tensors can be expressed in terms of the 



parallel propagator g aa '(x, x') between x and x' , and the 
two vectors n a = V a /i and n a > = V a '/i. Our two-point 
function is thus written as 

W aa '{x,x) = ai(a)g aa , + a 2 (u)n a n a > 

where, for spacelike separations (fi 2 > 0), ati and «2 are 
scalar functions of the geodesic distance n only. We will 
frequently use the following formulas [l[: 



9a n b' 

^a9bc' 

where 



A{g a b - n a n b ), 
~(A + C)(g ab n c 



9ab = g a 9c'b 

V a n b / = C{g a b' - 
+ 9ac>n b ) 



n a n b >) 



A = 



cot(/x/i?) 



and C 



csc(/x/i?) 



R ' R 
It is useful to introduce the quantity z given for spacelike 
separations (fi 2 >0or0<z<l)by 

-<»*(.&)■ <D 

All of our calculations are performed in the spacelike 
region < z < 1. They can be extended further by 
analytic continuation in the z variable to the cut com- 
plex plane C \ (l,oo). Finally, the Feynman propaga- 
tor G^ a , = i(Cl\TA a (x)A a i(x')\Cl) can be obtained as the 
limiting value of W aa > when approaching the branch cut 
z > 1 from above. 

We end our introduction by noting that several, quite 
interesting, recent works (see among many others) 
were devoted to the study of infrared effects in de Sit- 
ter. In these articles however, only massive (although 
interacting) fields are considered. IR pathologies are of 
course much stronger for massless fields and they appear 
already at the linear level. 

The R\ gauges. The gauge fixed action describing 
the massless vector field A a is 



Sx = J d^x^—g (\F ab F" b + \{V a A«) 2 ^ 



(2) 
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where F ab = W a A b - V b A a . 
motion is D ab A b = where 

(d-1) 



The resulting equation of 



Dab = 9abO - 



R 2 



g ab + (A - 1) V a V&. 



(3) 
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1. Flat space. In flat space (R —> oo), the Fourier 
transform of the Feynman propagator reads 



1 



1-A 



ie 



A (fc 2 - ie)- 



:k a k a ' • 



Using the (massless limit) of the scalar Feynman propa- 
gator (given, as always, for spacelike separations) 



G F (x,x r ) = 



d A k 



im Ki{mn) 



(2tt) 4 k 2 + m 2 - ie (2ir) 2 /i 



(4) 

(where K is the modified Bessel function of the second 
kind) and its derivatives with respect to x and x', we get 
the coordinate-space two-point function 



Woof = 



A + l 



8tt 2 A^ 2 



A- 1 



4tt 2 A^ 2 



n a n a > 



(5) 



For all values of A, this expression vanishes for large time- 
like and spacelike separations as expected. 

2. De Sitter space. Following [l|, we write the two- 
point function in de Sitter space as 



W£! (x, x') = a[ Ai (z) g aa , + < J (z) n a n a , . (6) 

where a\,a2 are scalar functions of the invariant quantity 
z defined in ([1]). The equation of motion D a bW h a , = 
implies two independent equations on ot\ and cx2- These 
equations are most easily solved in the Feynman gauge 
A = 1 and the result is found in 111 to be 



^\z) = 



1 



48ir 2 R 2 



1 



24ir 2 R 2 



1 - z 
1 

1 - - 

z 



- + i]ln(I 



i-i)ln(l-:) 



Instead of vanishing for largely separated points, the cor- 
relation goes to a nonzero constant function 



lim W { y\x,x') 



1 



24tt 2 R 2 



n a n a 



(7) 



This constant vanishes in the flat spacetime limit (i? — ► 
oo) as expected. The nonvanishing of the correlator in 
de Sitter space seems to contradict our experience with 
the more familiar flat space QFT. It is indeed a general 
result of the latter, sometimes referred to as the cluster 
decomposition principle Q , that the correlation functions 
must decay in the IR region. To our knowledge, only 
very special QFT models, e.g., confinement models like 
the massless scalar in two-dimensional flat space, seem 
to escape this result. 

As explained in the introduction, it turns out that this 
IR pathological behavior of the photon field is typical of 
de Sitter massless fields. The most notable examples are 
the mmc scalar and the graviton fields. The correlation 
functions of these fields exhibit even stronger pathologies 
as they grow in the IR region (see, for instance, [4J46|). 



This "IR problem" is at the heart of an open discus- 
sion in the community. It led some authors [6| to speak 
of a "striking violation of cluster decomposition proper- 
ties of the de Sitter invariant vacuum state." It has also 
been shown that this growing of the graviton twenpoint 
function renders scattering amplitudes divergent [9|. As 
a consequence, the authors of Q, among many others, 
suspect a quantum instability of de Sitter space, mean- 
ing that the de Sitter geometry is not a stable ground 
state of quantum gravity with a cosmological constant. 

Other workers in the field, mainly motivated by the 
fact that the correlations of some gauge invariant quan- 
tities fall off in the IR limit, claim that the nonvanishing 
of the two-point functions in the IR region is nothing but 
a gauge artifact. For instance, the electromagnetic field 
correlation function is well behaved in the IR region[l|: 



(F ab F 



a ' b ') - 1 1 (a [a ' a b ' ] +4r ir a [6 V^ 

} ~ 8n 2 R4(l-z) 2 3& ] +An ^] n J. 



The same is true for the graviton field as there are no 
growing terms in the correlation function of the Riemann 
tensor. 

In the following we will study in detail the IR problem 
in the case of the photon field in de Sitter. The IR be- 
havior of this field received considerably less attention in 
the literature than its counterpart for the mmc and the 
graviton field. This is certainly due to the importance 
of these two fields in inflation theory. It is also proba- 
bly due to the fact that in actual QED calculations the 
nonvanishing term ([7|) will not contribute to scattering 
processes (see the last section). This model will however 
prove to be precious in realizing that - at least for some 
situations - serious IR pathologies in de Sitter space are, 
without doubt, gauge artifacts. 

The equation of motion Q applied to the two-point 
function ((6]) leads after some algebra to a coupled system 
of differential equations. Hereafter we will only need one 
of these equations, namely, 



H 



hal-2 



A-A rf 
A rf -1 



' 2 " 

z 



1-A 



A - 1 



Oil 



-a 



1 - z 



a 2 



= 



(8) 



where the prime denotes derivation with respect to the 
d 3 

variable z, A<j = — — - and H is the hypergeometric op- 



erator defined by 



1 



H[a, 6, c]a = z{\ — z)a" + a' + ab a 

Another, moreover first order, relation between ot\ and 
a 2 is obtained by considering the two-point function of 
the field strength 



(F ab F a ' b ) 



a [ a ' £>'] i a W q'1 

A 9[a 9 b] + n [a g b \ n J 



3 



In fact, the equation of motion (in vacuum) V a F a b = 
implies the relation [l[ 



a 2 



8tt 2 R 2 (z- 1) 



+ 2(1 - z) (za[ + ai) (9) 



We now specialize our results to the d = 4 case for sim- 
plicity. The general solution in the Euclidean vacuum is 
found to be 



a[ X] (z) = 



1 



487r 2 i? 2 A 



1 



3(A+1) 3A-1 

— H h 

2(1 -z) 2z 

<^-l)(i + ^)Ml-,) 

3(A-1) 3A-1 



247r 2 i? 2 A 



(3A - 1) 



2(1 -z) 
1 



2-j 



2z 2z 2 



ln(l - z) 



and the IR behavior is given by 



lim W^,(x,x') 



z— >±oo 



247r 2 i? 2 A 



n a n a > 



(10) 



This result means that one can cure any IR bad behavior 
by going to the Landau gauge A — > oo. More importantly, 
we see that the IR behavior of the correlation function 
is pure gauge, thus asking whether or not the two-point 
function is well or ill behaved in the IR region is in itself 
-at least partially-a misleading question. 

Finally we note that there exists a special choice of the 
gauge-fixing parameter, namely A = | (more generally 
Ad in d dimensions), that cancels logarithmic terms and 
gives a particularly simple two-point function. We note, 
perhaps as a curiosity for the time being, that for large 
dimensions d, Xd goes to the Feynman gauge A = 1, which 
is the simplest gauge in fiat space. 

More pathological gauges. We now consider a 
larger family of gauge-fixing terms. Our motivation is 
that the IR pathology we exhibited in the last section, 
namely that the correlation function tends to a constant, 
might seem -I believe mistakenly- mild enough not to 
worry about. We will show that in this more general 
gauge, the two-point function is actually growing in the 
IR region, exactly like the mmc scalar or the graviton 
fields. We consider the action 



-g(S a A a )U{V a A a ) (11) 



plays this role: m 2 



■jp. This observation renders the 
introduction of a higher derivative gauge-fixing somehow 
more natural in the de Sitter case. Finally note that 
the use of a two-parameter family of gauge-fixing terms 
that are simultaneously nonzero does not overconstrain 
the system's dynamics. In fact, in the Faddeev-Popov 
procedure, the gauge-fixing term added to the action is 
merely the probability distribution P[ui] of the function 
to(x) that intervenes in the generalized Lorentz gauge 
condition \7 a A a = lj(x). This is of course a valid gauge 
condition and does not overconstrain the dynamics. 

The equation of motion reads 



ag ab + (A 



l)V a V 6 + -^DVaVft 
Tfl 



.4" 



(12) 



1. Flat space. The flat space Feynman propagator is 
found in Fourier space to be 



GL(fc) 



9aa' 

k 2 — ie 



+ 



1 - A 



A (fc 2 - ie) 2 
1 



i 2 A 



m 2 A 2 \k 2 — ie k 2 — — ie 



The reader will note the relative negative sign typical of 
higher derivative theories in front of the last propagator. 
Using (jU) and its derivatives the coordinate-space corre- 
lator is found in terms of Bessel functions. For all values 
of the gauge-fixing parameters it vanishes in the IR as 
expected. 

2. De Sitter space. We now consider the action (fTTj) 
in de Sitter space and thus set 1/m 2 = R 2 . Following the 
same methods described above we obtain the following 
equations on a\: 



„(4) 



4(t-2)(t + 3) 3(2z-1)(t 2 +t-16) , 



(z- l) 2 z 2 



-Oi\ 



12-(z-1)z(t-8)(t + 9) „ 
a 



(z-l) 2 z 2 1 
1 1 

_ V I 



(z- l) 2 z 2 " x ' \z z- 1 

(t 2 +t-2)(2(1-3A)z-3) 18-5r(r+l 



327r 2 Ai? 2 (z - 1)V 



327r 2 i? 2 (z- l) 4 



where we have added a higher derivative gauge-fixing 
term. Similar higher derivative gauge fixing terms are 
used in QCD, the electroweak theories, and were studied 
in 0. They were also used to study flat space pertur- 
bative quantum gravity 0]. We note that if we want A 
and £ to be dimensionless, such higher derivative gauge- 
fixing terms require the introduction of a mass quantity 
m . In flat space m is introduced by hand, while in de 
Sitter space the inverse of the de Sitter radius naturally 



where we introduced the parameter r 



Ob- 



taining this equation requires some elaborate algebraic 
manipulations. The calculations in this section were thus 
verified with the tensor algebra system xAct on Mathe- 
matica [loj . This equation can be solved in closed form. 
Using the asymptotic formulas of Legendre functions, the 
solution verifying (i) regularity at z — and (ii) flat space 
short-distance singularity is found. For noninteger values 



4 



of t it reads 



ax(z) 



2Q 2 (2z- 1) -7rcot(7rr)P r 2 (2z- 1) 
32ir 2 R 2 \ (z - l)z (r 2 + r - 2) 



3A - 1 
967r 2 i? 2 A 



(2z + l)ln(l-z) 



1 



32tt 2 j R 2 A(t- 1)(t + 2) 

1 (A + l)r(r + l) -2A 



(s-l)* 



Z - 1 



where P and Q are Legendre functions of the first and 
second kind, respectively. As before, a 2 is obtained 
by ©. Using the asymptotic behavior of the Legen- 
dre functions we find for instance that for r > 2 the 
two-point function grows like z T ~ l in the IR region. For 
(A = 0, £ 0) the two-point function can be expressed in 
terms of polylogarithms and the IR behavior is given by 



w aa , 



In z 



72n 2 £R 2 



We thus find that it grows logarithmically in the infrared 
region like the two-point function of the mmc scalar field. 
This IR growing is purely a gauge artifact since, as proven 
in the beginning of this Letter, no IR pathologies arise 
in the Landau gauge. That such a growing behavior is 
nothing but a gauge artifact for the photon field suggests 
that the situation might be the same for the mmc and 
the graviton fields. This possibility will be examined in 
detail in [nj. 

A physical decomposition. The previously de- 
scribed situation, namely, that the IR behavior of the 
two-point function in de Sitter space contains important 
gauge artifacts, makes it natural to look for a decom- 
position into physical and nonphysical parts. We show 
now that such a decomposition exists and is given by the 
following rewriting of the two-point function: 



(13) 



In fact, most of the physics is included in the "quantum'' 
action Wj defined by 



exp 



n 



= / £>L4 ]exp 



d i x y ^A a J a 



where J a is an external conserved current V a J a = 0. The 
tree-level expression is given by 

Wj = [dV x dV x , J a (x)W aa ,(x,x')J a '(x') 



where dV x = d 4 Xy/—g(x) is the invariant volume ele- 
ment. Since the current J a is conserved, integration by 
parts ensures that only /3i will contribute to Wj and it 
thus fully deserves to be referred to as the "physical" 
two-point function. An explicit form of the functions /3i 
and fa can be found by solving the system 

<*i=0i+ R 2 a 2 = (1 - z)& + z(l - z)ft. 

The solutions are unique if one requires regularity near 
z = and we obtain 



1 



16n 2 R 2 Vl 



1 



ln(l - z) 



9aa> (14) 



It is a remarkable fact that W^ ys is (A, £) independent 
and is well behaved in the IR region. In flat spacetimc, 
the Feynman gauge (A = 1,£ = 0) gives precisely this 
physical two-point function. The situation in de Sitter 
is quite different and no choice of the gauge-fixing 
parameters gives the physical two-point function. A 
natural question is whether a more general gauge-fixing 
term gives directly W^ ys (x, x') in de Sitter space. 
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